Quantum nondemolition measurements of a flux qubit coupled to a noisy detector 
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We theoretically study the measurement-induced dephasing caused by back action noise in quantum non- 
demolition measurements of a superconducting flux qubit which is coupled to a superconducting quantum in- 
terference device (SQUID). Our analytical results indicate that information on qubit flows from qubit to de- 
tector, while quantum fluctuations which may cause dephasing of the qubit also inject to qubit. Furthermore, 
the measurement probability is frequency dependent in a short time scale and has a close relationship with 
the measurement-induced dephasing. When the detuning between driven and bare resonator equals coupling 
strength, we will access the state of qubit more easily. In other words, we obtain the maximum measurement 
rate. Finally, we analyzed mixed effect caused by coupling between non-diagonal term and external variable. 
We found that the initial information of qubit is destroyed due to quantum tunneling between the qubit states. 

PACS numbers: 03.65.Ta, 03.67.Lx, 42.50.Lc, 85.25.Dq 



I. INTRODUCTION 

The concept of quantum nondemolition (QND) measure- 
ment which can be used for detection of weak force acting 
on system such as gravitational wave had been introduced to 
allow repeatable detection^i The basic requirement of a QND 
measurement is that a series of measurements of the variable 
of system should give predictable results.— It means that the 
observable A must commute with Hamiltonian H that de- 
scribes interacting system and detector, i.e., [A, H] = 0. A 
large number of works associated with QND measurements 
have been done in quantum opticr'— cavity QED^ and cir- 
cuit QEDii 

Recently, QND measurements have experimentally been 
achieved on superconducting flux qubit system by investi- 
gating the correlation between the results of two consecutive 
measurements]^ However, in general , the measuring appara- 
tus is a mesoscopic system (e.g., dc-SQUID) which is coupled 
to environment. Therefore, noise can affect discrimination of 
signal which reflects the information of qubit via coupling be- 
tween meter and qubit (see Fig. [TJb)). There exist a larger 
number of works in mesoscopic system where one is forced to 
think about the quantum mechanics of detection process^iiil^ 
and about the fundamental quantum limits which constrain the 
performance of the detectorJ^ii^ In practice, noise plays an 
important role in quantum measurement: quantum noise from 
the meter acts back on the system, such as qubit, at the same 
time, the information about the variable conjugate to the mea- 
sured variable is destroyed. This phenomenon is omnipresent, 
because information about system is carried away into the sur- 
rounding due to indirect environmental coupling via coupling 
to detector For a weak measurement, the detection may be 
quantum limited that the signal-to-noise of measurement, de- 
fined as the ratio of the amplitude of the oscillation line in the 
output spectrum to background noise, is no more than 4J^ 

In this article, we study the correlation between measure- 
ment probability and measurement-induced dephasing. Our 
model which is used to analyze the interaction between qubit 
and detector is based on the theory of QND measurement of 
flux qubit dynamics in a short time4^ The conclusion shows 



that the easier we gain information about one variable, the 
faster we lose information about its conjugate variable. It 
means that the dc-SQUID can be used as a quantum-limited 
mesoscopic detectorJi We also study the deviation of QND 
measurement due to the coupling between non-diagonal term 
of qubit and variable of detector. This article is organized as 
follows. In Sec. II the basic models were introduced firstly, 
from which we use master equation to discuss the dephasing 
caused by quantum noise and its correlation with measure- 
ment probability distribution. In Sec. Ill, we discuss the ef- 
fect of non-diagnose term on QND measurement. We show 
that the mix effect caused by coupling between non-diagnose 
term and detector can demolish the information of qubit. 



II. PROBABILITY DISTRIBUTION AND 
MEASUREMENT-INDUCED DEPHASING 

The Hamiltonian of system is given hy^ 

H{t) = Hs + Hrneter + Hint + Hdrive{t), (1) 

where Hg is the Hamiltonian of qubit which can be considered 
as a two level system 



rr e ^ ^ » 

Hs = -a-z + -^'^x, 



(2) 



where ct^ and are Pauli matrices, e and A are energy dif- 
ference and tunnel splitting of qubit. The Hamiltonian of de- 
tector which acts as an oscillator is 



(3) 



The Hamiltonian that describes coupling between qubit and 
oscillator is 



Hint = hga^a^a 



(4) 



from which we know that the variable that we want to measure 
is (Tz. g is coupling strength. In other words, we can conclude 
the state of qubit from the states of detector according to the 
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FIG. 1: (a) Schematic of flux qubit coupled to a SQUID which used 
for qubit readout as a detector, (b) Schematic of the process of in- 
formation detraction from qubit and noise injection from one-port 
detector 



value of (T^. We assign |0) to ~ 1 and |1) to cr^ = —1. The 
Hamiltonian of external driving of harmonic oscillator is 



Hdr^ve{t) = f{t)ia + a^), (5) 

where f{t) — 2/cos(ajrft) is the driving force with driving 
frequency uJd- For an ideal QND measurement, it is necessary 
that the condition [a^ ,H]^0 must be satisfied. However, 
for superconducting flux qubits, due to the off-diagonal term 
A 0, the measurement is not perfect QND. In case that 
A <C e is satisfied, ct^ still would be treated as a conserved 
quantity on the time scale determined by 1/e and one expects 
small deviations from an ideal QND case"". After neglecting 
the off-diagonal term and moving in the rotating frame with 
frequency ujd, we obtain the new Hamiltonian 

H ^ + h{duj + ga^)a^a + f{a + a^) (6) 

with doj — ujos — Wrf. Then we can get a QND measurement 
from the form of the Hamiltonian. In order to find the rela- 
tionship between measurement-induced dephasing and prob- 
ability distribution, we analyze the impact of noise on QND 
measurement in part A and make a more quantitative analysis 
in part B. 



A. THE IMPACT OF NOISE ON QND MEASUREMENT 

Based on the model mentioned by Ref. [16], we substitute 
a Gaussian white noise for the zero point fluctuation because 
of the fact that noise contributes a large part in the actual mea- 
surement. The distribution of qubit signal / conditioned on 
the qubit state = ±1 is 



1 



zexp 



2Siit ^ 



(7) 



where Xi{t) = -\/2Re[ai(t)], i = ±, Sn is noise spectra den- 
sity. After turning on the interaction with time t, the coherent 
state qubit-dependent amplitude is found to be^^ 



—iScJit — nt /2 



with qubit-dependent detuning 6uJi = 5uj- 
dependent amplitudes and phase given by 



(8) 

gaz and the qubit- 
(9) 
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FIG. 2: (color online) (a) Probability distribution considering zero 
point fluctuation for getting the measurement result for and the 
initial state |0)(0|, plotted as a function of detuning Suj. (b) Proba- 
bility distribution considering the back action noise instead of zero 
point fluctuation with the same parameter. The damping rate of res- 
onator K = 0. IGHz, the couple strength between qubit and oscillator 
g — 0.3GHz and the back action noise spectra density Sii = 2/k. 
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If we choose a time t 
becomesi^ 



1/e, then the signal difference 



da{t) w V2tA, 



(11) 



2i<po 



where A = /(e 
probability distribution 



^ ) / \/2. We obtain the outcome 



/crf(4^)(a.)o] 



(12) 



If we measure a rotated quadrature X^f, = (ae~"* +a^e*'^ ) / 
and choose the phase (f> — arg A. The qubit difference be- 
comes 5x{t) = \A\t, and the probability for the measurement 
is 



(13) 



In Fig.|2]we plot the probability of measuring the state 
under the condition that we consider the contribution of zero 
point fluctuation (Fig. |2](a)). We also plot the probability of 
measuring the state under the condition that we consider the 
contribution of back action noise instead of the zero point fluc- 
tuation during the measurement process (Fig. |2](b)). Due to 
the influence of back action noise, we can find that the strong 
measurement region, compared with the one that only consid- 
ered the zero point fluctuation, is reduced. It means that we 
need more time to reach the same measurement strength as the 
situation that we only consider contribution of the zero point 
fluctuation. 



B. CORRELATION BETWEEN MEASUREMENT 
PROBABILITY AND DEPHASING RATE 

We choose the initial state of system to be p(0) = po <E) 
i.e., at time t = the harmonic oscillator in the vac- 
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uum state. Then, we turn on qubit and oscillator interaction, 
after a period of time, qubit and oscillator entangle with each 
other At the same time, we perform a strong measurement 
of the flux quadrature x ^ a + a) . As a result, the pointer 
state (the oscillator state) is projected to the state If 
the measurement is strong enough that the states of oscillator 
\x±){x± I which correspond to the qubit states |0) and 1 1) re- 
spectively are orthogonal, the coherence of qubit is destroyed 
by the measurement. To describe this phenomenon, we must 
find out the equation of motion that can describes the system 
state evolution. The master equation used in circuit QED sys- 
tem is valid here3 The qubit density matrix obeys: 

p = Cp= --[H, p] + KV[a]p + iiV[a4p + fV[a,]p, 

where is damping superoperator defined by = 

{2LpL^ — L^Lp — pL^L)/2. In above expression, k, is de- 
cay rate of resonator, 71 is qubit relaxation rate and 72 is de- 
phasing rate. To figure out the effect induced by measure- 
ment due to coupling of qubit to the oscillator, we do not con- 
sider the coupling between qubit and environment. Therefore, 
we neglect relaxation due to coupled to environment and set 
7i = 0, (i.e., for an ideal QND measurement no quantum 
tunneling occurs between qubit states |0) and |1) during mea- 
surement). Eq. ( fT4] i can be solved under the condition that the 
density matrix with an expansion in positive P representation. 
In Ref. [10] the solution of the master equation was found to 
be 



aio(i) = aio(0)exp[-i(e - 172)* - i2g I a+{t)Q*_{t)dt'] 

Jo 

(15) 

and aoi(t) = a-io{t)- The amplitude of coherence state a±{t) 
are determined by 

a+{t) = al+cxp[-{K/2 + ig + iSuj)t]{a+{0)-al) (16) 



with a'"^ = —if /{iK/2 + ig + Soj), and 
a_{t) = al+exp[-(K/2-ig+i5w)t](a_(0)-al) (17) 

with ai = —if/{k/2 — ig + iSuj), where afj. is the steady 
coherence state value of the oscillator, and are am- 
plitudes of coherent state |q!-|_) and |q;_), which depends on 
qubit states, aoi is the coherence term. 

We can find that 72 is the intrinsic dephasing rate from 
Eq. ( fTsT i and r„i is the measurement-induced dephasing rate 
which is determined by for a project measurement and ex- 
pressed as followingiii 

n-)Kg^ 



r,„ = 2<7lm(«;ai) 



(18) 



where n-t 



\a%\^ = /V[kV4 + {6uj ± g)^] is aver- 
age number of photons in the resonator. We find that for 
Su; = ±\/g'^ + we will get a large probability if mea- 
sure the amplitude of the signal. In other words, when match- 
ing one of the frequencies of qubit (i.e., Slj = ig), the in- 
formation of qubit is encoded in the amplitude rather than 



phase which is the conjugate variable with amplitude. At the 
same time, the measurement-induced dephasing rate will get 
the maximum value. In Fig. [51 we plot the probability dis- 
tribution that measuring state conditioned on initial state 
is prepared on |0) (0| and measurement-induced dephasing as 
a function of detuning. When matching one of qubit fre- 
quency, both measurement-induced dephasing rate and mea- 
surement probability get a maximum value. In the region 
where crfdAjt) w 1, it corresponds to a strong projective 
measurement case. At the same time, the two coherent states 
of oscillator are well separated in phase space due to decay of 
coherence term at a measurement-induce dephaing rate r,„. 
In general case, at time t = 0, the system is in a product state 

l^(i = 0)) = i(|0) + |l))®|a), (19) 

where \a) is the initial state of oscillator At time t, the state 
of system can be written as 



|^(t))^-(|0)|a+) + |l)|a_)). 



(20) 



Then, if we know oscillator state is a;, the qubit state 
can be determined exactly. From Eq. ( fTSl l we find that the 
measurement-induced dephasing depends on the overlap of 
both two oscillator state aj^. In case that the oscillator states 
|a_) and |a+) are orthogonal, the measurement could be con- 
sidered as a strong projective measurement, i.e., in the region 
crf(|A|i) = 1. For the maximum dephasing value, the os- 
cillator states decay to a steady state significantly, which ex- 
plains the feature of Fig. |3]that both measurement-induced de- 
phasing rate and measurement probability get a high value for 
matching the qubit frequency. In this situation, the qubit state 
is encoded in the amplitude of coherence rather than phase. 
It means that the easier we gain information about amplitude, 
the faster we lose information about phase. The measurement- 
induced dephasing rate can also be expressed as 



|(a_(t)|a+(t))He 



(21) 



In addition, as shown in Fig. |3] with the value of damping 
rate k increasing, both structures of probability distribution 
and measurement-induced dephasing rate become flat with 
synchronous tendency, which consistent with the expectation. 
Increasing of driving force / only change the overall scale of 
(e)~(h) rather than the structure, which explains the trend that 
with the value of driving force increasing, the strong measure- 
ment region in (a)^(h) becomes wider and the measurement 
probability increases significantly, especially for matching the 
qubit frequency. 



III. DEVIATION OF A PERFECT QND MEASUREMENT 
DUE TO MEASUREMENT-INDUCED TUNNELING 

Taking into account the impact of non-diagonal term of 
qubit, we obtain the system Hamiltonian 



(22) 




Detuning So (GHz) 

FIG. 3: (color online) (a)~(d) Probability of measuring state for the initial state jO)(0| with damping rate of oscillator k = 0.1 GHz, 
K — 0.2 GHz, K = 0.3GHz, k = 0.4GHz, respectively, and qubit coupling strength g = 0.3GHz, t = 0.1ns, noise spectra density 
Sii = 2/k. (e)~(h) measurement-induced dephasing rate as a function of detuning Su between resonator and drive frequency with 
driving force / = IGHz. The parameters of (e)~(h) are same with those in (a)~(d), respectively. 



Then we can diagonalize the qubit Hamiltonian Hg. In new 
eigenvector space, the system Hamiltonian can be written as 



E 

meanwhile, the interact Hamiltonian can be written as 

Hint hguno^a, 



(23) 



(24) 



where E ^ v? + A^/2 is energy splitting of qubit and n 
is a vector that represents the direction of qubit basis before 
diagonalized relative to the energy basis 



(25) 



with rj = arctan (A/e). The eigenstates are denoted by the 
superposition of |0) and |1) 



|t)=cos||0)+sin||l) 



;) = -sin|| 



cos||l) 



(26) 



From Eq. ( l24l i we find that fluctuation induced by exter- 
nal noise in particle number h = a) a can causes transition 
between qubit states. Therefore, qubit acts as a spectrum ana- 
lyzer. We separate the interaction Hamiltonian into two parts, 
one is hg cos rjazCi^a and the other is hg sin riaxO^a. The first 
part does not affect the repeatability of QND measurement as 
discussed in previous sections. The second part will cause 
"spin flip" transition due to the coupling between qubit off- 
diagonal term and external variable. According to the equa- 
tion of motion, the time evolution of density matrix is obtained 
in interaction picture by 



in 



(27) 



We choose the initial state as the a product of qubit and 
detector, i.e., pifi) = /3o(0) ® Pd{0)- Then, tracing Eq. (|27] | 
over the detector degrees and expanding its right hand to the 
second order, we obtain the time evolution of reduced density 
matrix po- Eq. (|27] | becomes 



/^oW- -^y^ dT{[ai{t),ai{t-T)p\t-T)]{n'{T)n)n 



(28) 



After taking matrix elements between eigenstates of Hg, 
we finally obtain the equation of motion that can describe the 
reduced density matrix evolution 

Pik'it)^-f dTY,Mkk'wpiAt), (29) 
Jo 



in which we introduced the matrix 



Mkk'W — 

+ (nn^(r))[(a^(t-r)a^(i))z'..'4.-(a^W)/.-^(<^^(i-r))rfe']}. 

(30) 

To calculate the energy relaxation rate and dephasing rate, 
we compute the transition rate between qubit eigenstates by 
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computing evolution of the diagonal terms of reduced density 
matrix pn and poo, and the decay of the off-diagonal term 
poi with the initial state being a (Jx eigenstate. The transition 
rates between the two qubit eigenstates can be derived from 
Eq. 



Tj, = .g^sin^r/S'„„(a;oi) 
r-|. = SIT? 77S'„„(-woi), 



(31) 



where r-|- and are rates at which qubit is excited from 
ground state to excited state and decay from excited state 
to ground state, respectively, wqi = E/h, S'„„(w) = 
dre^"^(n(r)7i(0)). If the noise source is in the thermal 
equilibrium at temperature T, the transition rate must satisfy 
the balance condition F^/Fj^ = exp(— /3/icjoi) with effective 
temperature (3 = 1/kBTeff. It means that the qubit energy 
can be absorbed or emitted by detector and the ratio between 
positive and negative noise spectral density depends on the 
effective temperature. We can also obtain the expression of 
dephasing rate F^ 



where 9q = arctan (2(7/k) w S^/k is phase shift during 
measurement derived from Eq. ( fTOl i. which is identical with 
Ref. [10] and Ref. [19] with a coefficient of expression dif- 
ference that is caused by the difference of the lower limit of 
integral. It is quite a contrast to Eq. ( fTSl l that has double peaks 
for matching the frequencies of qubit. With the amplitude of 
off-diagonal term of density matrix increasing, the pure de- 
phasing rate become weak and relaxation will be dominant. 

Form Eq. ( |33] |, we find that the decay of coherence is caused 
by intrinsic dephasing and aggravated by measurement- 
induced dephasing. In addition, the deviation from the initial 
value of (cr^) become bigger as the time goes by due to tran- 
sition between the qubit states. Therefore, the state of qubit 
is acquired; meanwhile, the repeatability of measurement is 
destroyed. It means that it should not be considered as a QND 
measurement. Based on this effect, it can be used for "back 
action cooling" in many system, such as a cantilever coupled 
to a optical cavity^sS a noisy qubit coupled to a mechanical 
resonator,(author?)— nanomechanical resonator coupled to a 
driven superconducting resonantor— and cooling a micromir- 
ror by radiation pressure inside a optical cavity 



1 



r^ = -g sin ?y[S'„„(woi) + 5'„„(-woi)] 



9 • ? 

f sm ' 

+ cos^ rjSnni!^ = 0) 
_ Tt + F; 



(32) 



■70, 



where 7^ is qubit pure dephasing rate. For the case A = 0, 
the measurement is a perfect QND measurement which cannot 
causes transition between qubit states, in which the dephasing 
rate is 



dT{fi^ {T)n) 



D- 



(33) 



To calculate dephasing rate caused by fluctuation of particle 
number of oscillator during measurement, we should calcu- 
late correlator {6h{T)Sn), where 5n is quantum fluctuations 
around the mean of photons in the resonator For a driving 
resonator with frequency oj^, the following equation is valid''' 



a(t) 



d{t)], 



(34) 



where a is a classic part corresponding to the coherence state 
\a) and d{t) is a quantum part which can annihilate the coher- 
ence state to the vacuum state. Then the correlator becomes 



IV. CONCLUSION 

We have analyzed the QND measurement of a flux qubit 
coupled to a noisy detector The analytical results reveal 
that both of the measurement probability and measurement- 
induced dephasing have a similar trend that their values have 
peaks for matching one of the qubit frequencies and reduce 
at the resonance point, at which driving frequency equals the 
bare harmonic oscillator frequency. With the resonator de- 
cay rate increasing, both of the curves of dephasing rate and 
measurement probability become flat. Due to the coupling 
between the non-diagonal term of qubit and variable of de- 
tector, the noise of detector perturbs the initial state of qubit. 
As a result, the transition between the qubit states destroys 
the initial information of qubit. Therefore, there will be en- 
ergy exchange between qubit and detector according to the 
value of effective temperature of noise, which can be used 
for "back action" cooling. Conversely, the behavior of qubit 
which acts as a spectrum analyzer can also be used for ana- 
lyzing the property of noise. 



{Sn{T)Sn{0)) ^ (rf(r)d(O)) = ne' 



— A6ujt — i^t 



(35) 



where n is the average photo number in the resonator Under 
the condition that Slu = and weak coupling g ^ k, the 
measurement induced dephasing rate is 



Fm — KuOq, 



(36) 
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